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1. Introduction 

The nucleon mass provides an important benchmark calculation for lattice calculations. Fur- 
thermore, the chiral expansion of the nucleon mass is now theoretically well understood, having 
been worked out to two-loop order in the framework of SU (2) heavy baryon chiral perturbation 
theory [1]. This theoretical tractability and the relative cleanliness with which it can be calculated 
from the lattice, make the nucleon mass a candidate for scale setting. I would like to share with you 
some lessons I have recently learned about the nucleon mass, lattice QCD and chiral perturbation 
theory. In particular, by examining all the recent (some preliminary) dynamical lattice calculations 
of the nucleon mass, from the BMW [2], ETM [3], JLQCD [4], LHP [5], MILC [6], NPLQCD [7], 
PACS-CS [8], QCDSF/UKQCD [9] and RBC/UKQCD [10] Collaborations, I will show that chiral 
extrapolations of the nucleon mass have some previously hidden surprises. 

From the QCD trace anomaly, we know that the mass of all hadrons composed of valence light 
quarks, the up, down and strange flavors, are dominated by a quark mass independent term with 
chiral corrections which begin linear in the light quark masses. 



The exceptions to this rule are the eight pseudo-Nambu-Goldstone modes, ^ = {■K,K,r\}, whose 
masses would vanish in the chiral limit, and otherwise scale as 



The dots represent corrections to the masses which include terms polynomial in the quark masses 
as well as the much hunted for chiral logarithms and other chiral non-analytic functions. The 
unambiguous discovery of predicted chiral non-analytic behavior in hadron observables is held by 
many to be necessary and un-refutable evidence that our lattice calculations have entered the chiral 
regime, where the quark masses are light enough such that observables can be described by chiral 
perturbation theory (;tPT) [11] and systematically extrapolated to the physical quark mass limit. 

After reviewing heavy baryon X?T, I will examine the nucleon mass results from the above 
mentioned Collaborations and show they exhibit striking evidence of chiral non-analytic behavior. 
However, the observed quark mass dependence is not predicted by any theoretical understanding 
of QCD we currently have, as the nucleon mass results do not obviously scale as = + 
ntq + c'^i^n^q'^ + . . . as predicted, but rather as Mat = + a^i^y/m^, for pion masses in the range 
niji ~ 156-750 MeV. Given the variety of lattice actions used, as well as lattice spacings, lattice 
volumes and scale setting procedures, this behavior appeai^s to be a phenomenon of QCD and not a 
conspiracy of lattice artifacts. I conclude by exploring the precision which will be necessary to rule 
out this theoretically unmotivated quark mass dependence or accept it as a new mystery of QCD. 

2. Heavy Baryon Chiral Perturbation Theory 

I will begin with a review of heavy baryon chiral perturbation theory which was first developed 
by Jenkins and Manohar [12]. It is an effective field theory of non-relativistic baryons interacting 
with soft pions and photons, which provides a systematic expansion of baryon observables about 
the chiral limit. For simplicity I will focus the two flavor theory in the isospin limit. ' 

'For an earlier and more recent review of HB;(;PT, see Refs. [13]. 



Mh =Mi; +c'{mq + ... 



(1.1) 



m-^ =B{mqj +mq,) + ... 



(1.2) 
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When constructing an effective field theory of nucleons (and other baryons) difficulties arise 
immediately in developing a consistent power counting scheme. Including the leading interaction 
with the axial field, = d^K-T / (V^f) + the nucleon Lagrangian is 

^ = \j/N{i^-Mo)YN + gAWNYnY5A^¥N + --- (2.1) 
leading to the equation of motion 

i^YN = Mo\l/N, (2.2) 

In all higher dimensional operators in the effective theory, we are free to replace derivatives acting 
on the nucleon field with its mass. Phenomenologically, we know the chiral limit value of the 
nucleon mass is approximately the same size as the chiral symmetry breaking scale, Mq ~ A^,. 
So how does one consistently count powers of derivatives, {i^/K^)" ~ {Mq/Kj^)" ~ i^(l) when 
constructing the effective Lagrangian to a given order? 

This is the question answered by Jenkins and Manohar. When interacting with low energy 
pions and photons, the nucleon, in first approximation, can be treated as an infinitely heavy static 
field. One begins by introducing a new set of nucleon fields related to the original nucleons by a 
phase, similar to that used to convert the Klein-Gordon equation to the Schrodinger equation; 

A^v(x) = ^^>'^''^'->^(x). (2.3) 

Here, is the four- velocity of the nucleon field. In the rest frame, = (1,0) and the projector 
(1 +/) /2 ^ (1 + 7o) /2 becomes the non-relativistic projector. The Lagrangian then takes the form 

^=A^>.5A^,, + 2gAM<5-AA^,, + 2aMA^vA^vtr(^+) + ^(A^') +^(Mo ^) . (2.4) 

The mass term, \j/f^MoYN is now absent from the Lagrangian, replaced with an infinite tower of 
^{Mq"') suppressed operators, leading to a dual expansion. In Eq. (2.4), I have additionally in- 
cluded the leading quark mass dependent term with the spurion field 

J^+=^-(^^m^Q^ + t,'^mQt,'^y with ^2^I = exp(\/2/7r-f//) (2.5) 

and the quark mass matrix is mq = diag(wi,m) with 2m = m„ + The spin operator satisfies 
S^Nv = —3/4Nv and v • S = 0. With this phase redefinition, it is convenient to decompose the 
nucleon momentum as 

P^=Mov^+k^, (2.6) 

where is a residual soft momentum. Derivatives acting on the nucleon field now bring down 
powers of this soft momentum, d^Ny = ik^Ny, which vanishes for on-shell nucleons at rest. We see 
that the theory has a consistent power counting provided k^ ~ ^ . 

One can also explicitly include the delta degrees of freedom in the theory. However, this 
introduces a new mass parameter in the theory; the nucleon-delta mass splitting in the chiral limit 

Ao=Ma-M;vU^=o- (2.7) 

Phenomenologically, we know A ~ 290 MeV. In all present day lattice calculations, the pion mass 
is heavier than in nature, and so it is prudent to treat Aq ~ rtij^ in the power counting. One can 
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{a) (b) (c) 

Figure 1: The leading order (a) and next-to-leading order (b) and (c) contributions to the nucleon mass. 

package the delta fields in a totally symmetric flavor tensor and also utilize the Rarita-Schwinger 
representation for the spin. This leads to the Lagrangian^ 



^ = NivDN + 2aMNNtr{^+ ) - [iv • D - Aq] - 2YMT^' T^tr{ 

+ 2gANS-AN + 2gi^T^S-AT^+g^{T^A^N+NAn^) . (2.8) 

The axial couplings, gA, gm and ^aa are phenomenologically very important and these axial in- 
teractions give rise to the leading non-analytic quark mass dependence of the baryon masses. For 
example, the nucleon mass to leading loop order, depicted in Figure 1 is 



Mn=Mo- 2aM(At)< - TJ-^K - ^zr^=^("i;r, Ao, Ai) (2.9) 



where I have defined the standard chiral log function as 



^(«z,A,M) = iA'^m^ + ie)yHn ( ^±^^^^] _ 1^2,^ ('^) -A^ln('4] ■ (2.10) 

The first term in Eq. (2.9) is the chiral limit value of the nucleon mass while the second term, 
denoted as the LO nucleon mass correction, comes from the tree graph. Fig. 1(a). The third and 
fourth terms are the NLO corrections to the nucleon mass arising from the one loop graphs depicted 
in Fig. 1 (b) and (c) respectively. Additionally, the explicit inclusion of the delta degrees of freedom 
leads to an infinite renormalization of the LEC aM, as well as a finite renormalization of and 
Mq. Furthermore, because the mass parameter Aq is a chiral singlet, it leads to an a priori unknown 
finite renormalization of all LECs in the theory [14], for example 

Mo^Mo + ^^+<^, aM(M) = «M-2gLv^lnf4V^^+<'''' 



gA^gA(l+<VA;t + ...), gAN^gAN{^+df^Ao/Ax + ...) . (2.11) 

As we do not have the ability to dial the nucleon-delta mass splitting Aq, as we can the quark 
masses, it is prudent to not keep explicit track of this finite renormalization but rather absorb it in a 
redefinition of the LECs. In this way one accounts for all finite polynomial dependence upon Aq. 

The price one pays for this formalism is that Lorentz invariance is not manifestly preserved. 
However, the symmetry is straightforward to recover order-by-order in inverse powers of Mq. The 
method is known as reparameterization invariance [15].^ Consider a small shift in the the nucleon 
momentum, Eq. (2.6), such that 

V — >v + e/Mo, k — >k-e. (2.12) 



^The velocity subscript on the baryon fields is implicit here and following. 
^For an equivalent formulation see Ref. [16]. 
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The nucleon momentum is invariant under this shift. Requiring the S -matrix elements of the theory 
to also maintain invariance under this reparameterization recovers the Lorentz invariance order-by- 
order in 1 /Mq. One manifestation of this requirement is a constraint on the coefficients of certain 
operators in the effective Lagrangian. Specifically, the coefficients of certain higher dimensional 
operators are exactly determined from a connection with lower dimensional operators. For exam- 
ple, considering the kinetic operators for the nucleon, one finds 



^ = Niv-dN- 



'N,^N + cN \,/ N+> 



2Mo 



2M( 



(2.13) 







The coefficient of the first recoil operator is exactly constrained to be — 1 while the coefficient of 
the second recoil operator (l/2Mo)A'^(v ■d)^N is unconstrained. However, there is still a freedom 
to make a field-redefinition to convert the Lagrangian into the following form [17], 



^ = NivdN-N 



2Mo 



N 



with 



(2.14) 



This choice provides the familiar- form of the non-relativistic propagator, provided the reciol oper- 
ator is re-summed to all orders. 



v-k + ie 



i{l+f)/2 



vk- 



2Mo 



(2.15) 



Calculating the nucleon mass to next-to-next-to-leading order (NNLO) requires this new recoil 
operator as well as other fixed-coefficient operators related to the axial coupling and a host of 
operators with unconstrained coefficients. The calculation is still a one-loop calculation, but with 
operator insertions on the internal nucleon (delta) lines. With the explicit inclusion of the delta 
degrees of freedom, the nucleon mass is given at NNLO by [14] 



Mn= Mo-2aMil^)ml- 



In 



(47rA)2'"'^ 

3fcA(M) 



6aM(M) 



'i 81 + 58' 



AN 



?In 



2Mo 



8gAw(«M(M)-yM(M)) 2 

(47r/,)2 



{4nf^)2 4Mo(47rA)2 32Mo(47rA)2 
where the new chiral log function is defined as 



m^,Ao,IJ.) 
(2.16) 



ie In 



A+VA2 



A- VA2 



m^ln 



le 



2A2ln 



4A^ 



(2.17) 



The formula for the nucleon mass is in fact known to next-to-next-to-next-to-leading order (NNNLO) 
in a theory without explicit deltas [1]. The heavy baryon theory has also been extended to quenched 
XVT [18] as well as partially quenched [19] in which the baryon masses have been determined 
to NNLO [14]. In the framework of heavy baryon ;tPT, the leading and sub-leading lattice spac- 
ing corrections to the baryon masses have also been determined for Wilson fermions [20] , twisted 
mass fermions [21], domain- wall on Wilson [22], domain- wall on staggered [23], domain-wall on 
anything [24] and staggered on staggered fermions [25]. 
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P P-k P-q 

{a) {b) 

Figure 2: Perturbative graphs contributing to e+e^ scattering, (a). A one-loop contribution to the nucleon 
scalar form factor, (b). 

2.1 Scalar Form Factor of the Nucleon 

The construction of heavy baryon x^T^ begins with the nucleon as infinitely heavy static fields. 
Therefore, a naive use of the theory will lead to incorrect analytic structure in the formula for vari- 
ous observables. Static quantities such as the nucleon mass, will be insensitive to these problems in 
the range of convergence of the theory. However, dynamic quantities with external momentum in- 
sertions, such as pion-nucleon scattering, the scalar form factor, etc. need to be analyzed with more 
care. In some cases, for kinematic reasons, the power counting of the theory must be rearranged. 
This is not a new phenomenon. To provide a familiar example, consider e^e^ scattering near the 
Z-pole. At any finite order in the perturbative calculation of the cross section, see Fig. 2{a), the 
sum will fail to reproduce the observed cross section. Re-summing the self-energy of the Z will 
generated a width, leading to the correct form for the cross section. 



For a specific example relevant to heavy baryon ;^PT, I will consider the scalar form factor of 
the nucleon, defined by the matrix element,^ 

a{t=q^) = {N{P)\mqqq\N{P-q)), (2.19) 

In Fig. 2{b), we depict one of the two one-loop diagrams contributing to this matrix element, the 
other having an internal delta state. It is well known that the evaluation of this matrix element using 
the static nucleon propagator leads to an unphysical singularity. So with hindsight in hand, I lets 
begin by re-summin the leading recoil coiTcctions to the nucleon propagator, Eq. (2.15). One then 
finds this loop contribution to the form factor 



where h{S, k) is a function of the spin operators and momenta. The exact form of h is not important 
for the manipulations I will make, as they will only modify the renormalization prescription. First 
recall the coefficient of the operator (l/2Mo)A^(v • d^N is not constrained as we are free to make a 
field redefinition to change it. Therefore, let us make the change P — > in the nucleon propagator. 



"^The scalar form factor was first computed in heavy baryon ;|jPT in Ref. [26] in tlie SU (3) theory. 
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I will also factor out the 2Mq from the nucleon propagator, and with some malice aforethought, I 
will make a slight adjustment to h{S,k), leaving us with 

^nIo(2. Sl f [Mof-l^+Mo]h'{S,k) 

P .1 {2nY {2MQvk-B){{k-qY-ml){B-ml)- ^ ' 

For on-shell nucleons, Mof = ^, such that the form factor takes the form 

^nlo,2. 8a f ^ It-HMo h\S,k) 

'^'^ pj {2KYMl-{P-kY{{k-q)^-ml){B-ml)- ^ 

Up to the numerator structure h'{S,k), this is exactly what we would have gotten by starting with 
the original relativistic Lagrangian, Eq. (2.1). Becher and Leutwyler, in 1999, wrote a nice paper 
related to these issues [27]. In addition to formally demonstrating the equivalence I have just dis- 
cussed, they further showed how to regulate the integrals appearing with the re-summed nucleon 
propagator, as in Eq. (2.22), in a manner which separates the chiral non-analytic pieces (the in- 
frai^ed pieces) from the polynomial in m^y pieces (the regular pieces). This has come to be known 
as covariant baryon xPT. Furthermore, the kinetic recoil corrections of the nucleon remove the 
unphysical singularity which occurs with a naive application of heavy baryon ;^PT, as I shall now 
demonstrate. Evaluating the integral in Eq. (2.22) using the rules of Ref. [27], one finds the scalar 
form factor is given at NLO by 




in 1+ "^^Z^'^' 



\/l-^/(4'4 




. (2.23) 



Near the pion-production threshold, t = 4m^, there are singularities in both terms of this NLO 
contribution to o{t), which exactly cancel. In a naive application of heavy baryon %VY, one would 
not have kept the second term, counting m-nj (2Mo) ^ 1 moving this term to NNLO and resulting in 
an unphysical singularity at ? = 4/m| from the first term alone. But as we have seen here, near this 
kinematic threshold, the power counting of the theory must be re-arranged to maintain a consistent 
and correct theory. Near t = 4wi|, one should count 



2Mq " V ' 4wi^ 



1 - ^ . (2.24) 



Working away from the kinematic thresholds such as this, one is using the power counting 

mn ~ Mo . (2.25) 

Furthermore, working away from these kinematic thresholds and not explicitly including the delta 
degrees of freedom is equivalent to choosing the power counting 

mn ^ Mq <S^Ma-Mn. (2.26) 

In summary, covariant baryon ;^PT is not a new effective field theory, but is heavy baryon ;t;PT 
with a re-summed class of diagrams [27]. The choice to use the original formulation of Jenkins 
and Manohar or to use the new covariant formalism of Becher and Leutwyler is a matter of taste. 
As with all theories, the use of heavy baryon ;fPT, or the so derived chiral extrapolation formulae, 
requires consumer cai^e; the naive use of either formalism may lead to inconsistencies. 
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NNLO - ml, with gA=l-2(l),giN=1.5(3) NNLO - mj, with gA=1.2(l),giN=0 




0.8 I ' ' ' ' ' ' ' ' ' ' ' ' ' ' ' 1 I ' ' ' ' ' ' ' ' ' ' ' ' ' ' ' 1 

0.0 0.2 0.4 0.6 0.8 '.0 0.2 0.4 0.6 0.8 



gA=1.2(l), gAN=l .5(3) gA=l .2(1), gAN=0 




0.0 0.2 04 0.6 0.8 0.0 0.2 0.4 0.6 0.8 



mjil^nk) (^) m, 1(242 nk) 

Figure 3: Chiral extrapolation of LHP C nucleon mass results using the NNLO continuum formulae with 
(a) and without (b) explicit delta degrees of freedom. In (c) and {d) we plot the corresponding contributions 
to the nucleon mass order by order 

3. Comparing with Lattice QCD Results 

With a review of heavy baryon ;^PT in hand, I will now turn to comparing this formalism to all 
the recent dynamical lattice calculations of the nucleon mass. I will begin with a review of the LHP 
results [5], with which I am the most familiar. The LHP Collaboration performed a mixed-action 
calculation^ with domain-wall valence fermions on the coarse, a ~ 0.125 fm, L ~ 2.5 fm Asqtad 
improved rooted staggered gauge ensembles with pion masses in the range niji ~ 290-750 MeV. 
The relevant mixed-action effective theory for this calculation is well established [23, 24, 29]. 
However, there were too many unknown pai^ameters in the mixed action formula to perform a 
successful chiral extrapolation. Therefore, continuum SU (2) heavy baryon ;t;PT, including explicit 
delta degrees of freedom was used to perform the chiral extrapolations resulting in 

Mw = 954 ±42 ±20 MeV. (3.1) 

The first error is statistical while the second error is a systematic error resulting from varying 
the LECs which were fixed in the extrapolation, for specific details on the extrapolation analysis, 
see Ref. [5]. The resulting chiral extrapolation is depicted in Fig. 3(a), plotted as a function of 
the pion mass over the approximate chiral symmetry breaking scale, A.^ ~ l^/lnfa where /o = 
121.9 MeV [30]. A similar extrapolation was performed without deltas (gA/v = 0), resulting in the 
extrapolation depicted in Fig. 3(^7), which is qualitatively the same. In Figs. 3(c) and {d) I plot the 

^The strange quark propagators and the majority of the Ught quark propagators for this calculation were performed 
by the NPLQCD Collaboration [28], 
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NNLO Covariant, with = 1 .2( 1 ) Ma, = o-^ + m^, 




0.0 0.2 0.4 0.6 0.8 o.O 0.2 0.4 0.6 0.8 

Figure 4: Nucleon mass extrapolation with full one-loop covariant formula, (a). A theoretically unmoti- 
vated extrapolation ofMj^ linear in (b). 

resulting order-by-order contribution to the nucleon mass from these two fits respectively. In these 
plots, the arrows denote the values of the pion mass at which the calculation was performed. 

For comparison purposes, in this talk, I have also performed a chiral extrapolation using the 
full one-loop covariant baryon ;t;PT formula, which can be found in Ref. [31]. This extrapolation is 
depicted in Fig. 4(a), which is again found to be qualitatively similar to the heavy baryon ;^PT fit. 
The 68% error band on the covariant fit is smaller than on the heavy baryon fits. This is the result 
of more LECs being constrained, and not the result of the quality of the fit. 

Despite the questionable convergence of heavy baryon ;^PT for these masses, the theory is 
consistent both with the lattice results and the physical point. However, these fits all required phe- 
nomenological input. Most notably, the axial couplings have to be fixed in these extrapolations, 
otherwise when left as free parameters, gA and g/^ float significantly away from their phenomeno- 
logically known values. Furthermore, it is impossible not to notice the incredible linearity in the 
nucleon mass results as a function of niji. In fact, performing an extrapolation with a theoretically 
unmotivated linear fitting function (depicted in Fig. 4{b)), one finds 

MA, = a^ + afm;r = 938±9MeV! (with ;t:2/dof= 1.46) (3.2) 

To highlight how wrong this fit ansatz is, consider the pion-nucleon sigma term, which is the scalar 
form factor discussed in Sec. 2.1, at zero momentum transfer, 

d 

OnN = a(0 = m,^3 — Mn . (3.3) 
aniq 

Chiral symmetry dictates that OjiN/niq const, in the chiral limit, while for this fit function, 
Eq. (3.2), O-jiNhnq — > oo. Still, one can not ignore the appai^ent phenomenological success of this 
fit in describing the LHP lattice results. Is this a bizarre conspiracy of lattice artifacts or a physical 
phenomenon? As the LHP calculation was performed at a single lattice spacing, and a single 
volume, one must examine the numerical calculations of other lattice groups to study this question. 

3.1 Comparing with A^y = 2 + 1 calculations 

I will begin with an examination of the results from A'y = 2 + 1 calculations, those of the BMW 
Collaboration [2], the MILC Collaboration [6], the NPLQCD [7], the PACS-CS Collaboration [8] 
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Table 1: Details of the Nf = 2+1 dynamical lattice calculations. 



Collaboration 
RBC/UKQCD 
PACS-CS 
NPLQCD 
MILC 
LHP 
BMW 



a[fm] 
0.114 
0.091 
0.125, 0.09 
0.125, 0.09, 0.06 
0.125 
0.125,0.085,0.065 



L[fm] 
2J4 
2.90 
2.5 

2.5, 3.0, 3.5 
2.5 

2.0, 3.0, 4.0 



m-nL > 
4.56 
2.29 
3.7 
3.7 
3.7 
4 



scale setting 

Ma- 

Ms, Mfl- . . . 



and the RBC/UKQCD Collaborations [10]. In Table 1 1 list relevant parameters from these lattices. 
As can be seen, there is a wide range of pai^ameters used as well as several independent methods of 
scale setting. In Figure 5, 1 plot the resulting values of the nucleon mass from these collaborations, 
in all cases using the scale setting method of the given Collaborations to convert all results to GeV. 
When it does not crowd the plot too much, I also add the LHP results as a reference. 

In Fig. 5{a), one observes the coai^se (a ~ 0.125 fm) MILC results are systematically higher 
than the LHP results. This is most likely an additive lattice spacing correction to the nucleon 
mass. The super-fine (a ~ 0.06 fm) MILC results are very consistent with the LHP results with 
the lower error band coming from a fit to the LHP results alone. The super-fine MILC results, 
which are preliminary, are consistent with this linear in niji extrapolation down to niji 220 MeV. 
In Fig. 5{a), I have also added a preliminary result from NPLQCD on the fine {a ~ 0.09 fm) 
MILC lattices with domain-wall valence fermions, which is also consistent with the LHP results. 
This is suggestive that the discretization errors in the mixed-action domain-wall valence on MILC 
sea fermions are small. In Fig. 5(^) I have plotted the preliminary RBC/UKQCD results which 



MILC: MN = a^ + Q-f ni„ 




(a) 



> 
O 



4 NPLQCD 
MILC: coarse 
MILC: super fine 
LHP 



n,/(2^/2 7:fo) 
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■ m 








1.2 
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i PACS-CS 






i LHP 



^/ i RBC 
/ i LHP 




(b) 


0.2 0.4 0.6 

mJ(2^/2J:fo) 













i i BMW: a~ 

It 
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-0.125 fm 
'0.085 fm 
'0.065 fm 
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0.2 0.4 0.6 

m,/(2V2";r/o) 






Figure 5: Nucleon mass results from all Nf = 2+1 flavor lattice calculations. 
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are an update of those in Ref. [10]. One observes the Ughter two points are consistent with the 
LHP results while the heavier two points are inconsistent as measured with the statistical errors. 
However, RBC/UKQCD use an independent method of scale setting, which may account for this 
discrepancy. In Fig. 5(c), 1 have plotted the recent results from PACS-CS [8], which are also 
consistent with those of LHP. The eiTor band in Fig. 5(c) is the result of fitting the lightest four 
PACS-CS results. Although their lattice spacing is reasonably small, the lightest PACS-CS result, 
at the impressive niji ^ 156 MeV has a rather small value of m^L and is thus expected to have a 
significant volume correction, although perhaps not bigger than the current statistical error bar. The 
BMW Collaboration has the most impressive set of results, plotted in Fig. 5 {d) . They have results 
at three different lattice spacings, and all points plotted have m^jL > 4, with the lightest pion mass 
m„ ~ 190 MeV. Within the statistical errors, there are no discernible discretization errors. 

Given the large number of independent results from the various collaborations, and with the 
three independent scale setting methods, we can take the trends of the nucleon mass results as a 
function of mji seriously. The most striking feature is that all these Nf = 2-\- \ results display 
this remarkably linear dependence upon the pion mass. For the range of pion masses at which 
these calculations were performed, it would seem that QCD has conspired to produce this highly 
unexpected chiral non-analytic behavior of the nucleon mass, such that Mj^ ~ + cCy ,yjn^. Before 
discussing what is needed to rule out this theoretically unmotivated chiral behavior in favor of 
that expected from heavy baryon ;tPT, I will examine the results of the Nf = 2 flavor dynamical 
calculations. 

3.2 Comparing with Nf = 2 calculations 

1 now make a similar examination of the results from the Nf = 2 calculations from the ETM [3], 
JLQCD [4] and QCDSF/UKQCD [9] Collaborations. In Table 2, 1 collect pertinent parameters 
from the different lattices used in the calculations. In Fig. 6, 1 plot the resulting values of the 
nucleon mass, using the scale setting methods employed by the respective groups. The Nf = 2 
calculations all have smaller volumes (and minimum values of m^jL) than the Nf = 2+1 calcula- 
tions. There is also less independence in scale setting methods; all three groups use ro [32] while 
ETM also uses fj^ to set the scale. As compared to the A'y = 2 + 1 calculations, the Nf = 2 calcu- 
lations are more varied with respect both to each other as well as with respect to the coarse LHP 
and super-fine MILC results which 1 am using as a reference. With either scale setting method, 
the ETM results (continuum extrapolated), depicted in Fig. 6{a), are. systematically higher than the 
LHP/MILC results, with the results using fj^ to set the scale [33] being the most different. The 
ETM results have been computed in smaller volumes (at smaller mj^L) and the predicted one-loop 
volume corrections [34] will lower the results at the few percent level, bringing them in closer 
agreement with the LHP/MILC results. The resuhs of the JLQCD Collaboration (a ~ 0.118 fm 



Table 2: Details of the Nf = 2 dynamical lattice calculations. 



Collaboration 


a[fm] 


L[fm] 


my[L > 


scale setting 


ETM 


0.083, 0.0655 


2.10, 2.66 


3.3 




JLQCD 


0.118 


1.9 


2.8 


ro 


QCDSF/UKQCD 


0.085-0.067 


1.5-2.6 


2.1 
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Figure 6: Nucleon mass results from all Nf — 2 flavor lattice calculations. 

clover fermions) are depicted in Fig. 6{b). The raw numbers are in approximate agreement with 
the LHP/MILC results. However, the JLQCD results were calculated with smallish volumes. They 
adjusted the results for predicted finite volume corrections, whcih I have listed as JLQCD: FV in 
Fig. 6{b). One observes these numbers at the lightest two points are not in agreement with the 
LHP/MILC results. Furthermore, they display definite curvature (non linearity in mji behavior). 
The QCDSF/UKQCD results {a ~ 0.085-0.067 fm) are expected to have the largest finite volume 
corrections, for which these numbers have not been adjusted. The results are from Ref. [9] as well 
as preliminary new numbers. The (new) preliminary numbers from QCDSF/UKQCD are not in 
agreement with the LHP/MILC results, and the lighter points are expected to have larger finite vol- 
ume corrections, possibly as large as 5-10%. However, the results from Ref. [9] are in agreement 
with the LHP/MILC numbers. 

Although the Nf = 2 nucleon mass results are in less agreement with the LHP/MILC results 
(as compared to all the Nf = 2+1 results which are all in good agreement), aside from the JLQCD 
finite volume adjusted results, the Nf = 2 results are also consistent with the theoretically un- 
motivated linear in behavior. There appears to be overwhelming "lattice phenomenological" 
evidence that this unexpected chiral non-analytic behavior is a real QCD phenomenon and not sim- 
ply a bizarre coincidence of lattice artifacts. A natural question to ask is what level of precision 
is required in our numerical nucleon mass results to rule out this possibility? I will address this 
question in the next section. 

4. Challenges with chiral extrapolations for baryon observables 

Given the numerical state of affairs, in particular the nucleon mass results from the Nf = 2+1 
calculations, I can only make qualitative statements about the needed precision to rule out the 
straight-line fit in favor of the heavy baryon ;^PT analysis. To quantify my qualitative statements, I 
will use a /dof > 2 to rule out a particular fit. I will continue to use the LHP results [5] as my 
reference. I will perform the following excercise; by hand I will lower the lightest mass results by 
3%, raise the third lightest point by 2%, lower the fifth lightest point by 1% and leave the second, 
fourth and heaviest points as they ai^e. In this way I am adding curvature to the results which can be 
accomodated with the unknown LECs in the heavy baiyon ;^PT extrapolation formula, but which 
will rule out the linear in niji fit. See Fig. 7(a). One can then ask how small must the error bars 
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Figure 7: Curvature added to the LHP results, (a). £TM, and LHP nucleon mass results in ro units, {b). 



be on the nucleon mass results to rule out the straight-line in m-n fill The answer is 1 % which 
returns a 'X^ jdof = 2.08. With 2% error bars, the jdof = 0.52, in both cases with four degrees 
of freedom. This means we would have to know each nucleon mass point at the 1 % level including 
all sources of systematic error! A challenging endeavor to say the least. 

At this level of precision, all systematic errors may come into play; the lattice spacing correc- 
tions (present in all but perfect lattice actions), finite volume corrections, heavy quark corrections 
(from not having the c, h and t dynamical quarks), scale setting errors, etc. From the analysis in this 
talk, it appears that scale setting may introduce one of the larger uncertainties. All the A'/ = 2 + 1 
calculations are in good agreement with each other, and have used three independent scale setting 
methods. However, the ETM results have been determined with two independent scale setting 
methods, the use of ro and also the use of f-^. From Fig. 6(Zj), one can see a striking difference 
between the two methods which should in principle agree. If one instead plots the LHP and ETM 
results in ro units, the disagreement between the two sets of results becomes even less pronounced, 
although still with a slight systematic shift. See Fig. lib). Given the lightest quark mass results 
of the BMW and MILC Collaborations, it seems evident that the resolution of this puzzle does 
not lie with lighter quark masses, at least those which are light but still heavier than the physical 
pion mass. Collectively, we will need to address all the known lattice artifacts with numerical cal- 
culations at finer lattice spacings, larger volumes, and perhaps more care in addressing the eiTors 
associated with scale setting and comparing scale setting methods. 



5. Conclusions 



The nucleon mass has often been thought of as an important and clean benchmark observable 
with which to calibrate ones lattice calculations, and a candidate quantity for scale setting. By 
clean, 1 mean relatively free of systematic uncertainties. Little did we know that (lattice) QCD 
has conspired to make the nucleon mass fai^ more mysterious and wrought with uncertainty than 
previously thought. For years, we as a lattice community have been searching for the tell-tail 
signs of the chiral logs (and other predicted chiral non-analytic functions). A close examination 
of calculations of the nucleon mass reveals that the results from all dynamical A^y = 2 + 1 and 
some of the 'Nf = 2 lattice calculations can be very well described by the chiral non-analytic (but 
unexpected) function 

M/v = + aJ^niTt ■ 
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Perhaps more shocking, extrapolations of the lattice results to the physical pion mass produce a 
number in striking agreement with the physical nucleon mass; from the LHP results one finds 
Mj^ = 938 ±9 MeV. Given the large amount of independent lattice actions, lattice spacings, lattice 
volumes and scale setting methods, this appears not to be a bizarre conspiracy of lattice artifacts 
but rather a phenomenon from QCD. Resolution of this puzzle (ruling out this linear in rrij^ function 
in favor of the heavy baryon ;^PT formula) may require nucleon mass results with uncertainties at 
the sub 1% level, including all systematic errors. 

One thing is now clear, the nucleon mass is a bad observable to use to hunt for predicted 
chiral non-analytic behavior. A determination of the phenomenologically interesting [35] pion- 
nucleon sigma term, Ojin, I predict will further require lattice calculational results on both sides 
of the physical point. From the LHP numbers, one finds an extrapolation including explicit delta 
degrees of freedom results in a^N ~ 84 it 17 it 20 MeV. The result without explicit deltas returns 
a value closer to the current consensus, o^ij^ 42 it 14 it 9 MeV. Using the linear in mj^ fit, which 
seems further justified by the BMW and preliminary super-fine {a ~ 0.06 fm) MILC results, returns 
a value of o^n ~ 67 it 2 MeV, which is in good agreement with the determination of the ETM 
Collaboration [3]. 

As a final note, I want to mention that all the heavy baryon ;^PT analyses of nucleon mass 
results have required one to input the value of the axial coupling, generally with the value ~ 1 -2. 
Although possibly acceptable, this is distasteful. One would really like to see the value of this 
coupling emerge from a fit to the nucleon mass itself, which would be a smoking gun in a hunt for 
predicted chiral non-analytic dependence upon the light quark masses. However, given the state 
of the lattice results for the nucleon mass, see Fig. 5, this seems quite unlikely to happen. As an 
alternative, one should move in the direction of global fits, which include both the the nucleon 
mass and the nucleon axial charge, as well as other observables which can help to constrain the 
phenomenologically more interesting couplings in the nucleon Lagrangian. 
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